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Abstract. We prove the spacetime positive mass theorem in dimensions less than 
eight. This theorem states that for any asymptotically flat initial data set satis- 
fying the dominant energy condition, the inequality E > \P\ holds, where {E,P) 
is the ADM energy- momentum vector of the initial data set. Previously, this the- 
orem was proven only for spin manifolds [331 . Our proof is a modification of the 
minimal hypersurface technique that was used by the last author and S.-T. Yau 
to establish the time-symmetric case [551152]. Instead of minimal hypersurfaces, 
we use marginally outer trapped hypersurfaces (MOTS) whose existence is guar- 
anteed by earlier work of the first author [12] . Since MOTS do not arise from a 
variational principle, an important part of our proof to introduce an appropriate 
substitute for the area functional used in the time-symmetric case. As part of our 
proof, we establish a density theorem of independent interest that allows us to re- 
duce the general case of our theorem to the case of initial data that has harmonic 
asymptotics and satisfies the strict dominant energy condition. A refinement of 
the density argument allows us to approximate any given data set by one which is 
identical with a Kerr slice outside a compact set preserving the dominant energy 
condition. This enables us to give an alternative proof of the main theorem by 
reducing it to the positive energy theorem. 



1. Introduction 

The main result of this article is the following theorem. 

Theorem 1 (Spacetime positive mass theorem). Let 3 < n < 8, and let {M,g,k) be 
an n- dimensional asymptotically flat initial data set satisfying the dominant energy 
condition. Then 

E > |P|, 

where {E,P) is the ADM energy-momentum vector of {M,g, k). 
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Precise definitions of all technical terms used in this article are given in Section [2j 
We briefly survey earlier results: A particularly important special case of this 
theorem is when k = 0, which is called the time-symmetric case, or sometimes the 
Riemannian case. In this case, we have P = 0, and the dominant energy condition 
reduces to the assumption of nonnegative scalar curvature of g. The last author 
and S.-T. Yau proved the time-symmetric case of the theorem in dimension three in 
two articles from 1979 and 1981 [25l 126] . In [2H] they generalized their proof of the 
time-symmetric case to dimensions less than 8 without any topological assumption, 
as explained in detail in [22]. In 1981, they considered the general case /c ^ in 
dimension three and succeeded in proving that E > Ohj solving Jang's equation [27j. 
Later, E. Witten discovered a completely different proof that E > \P\ in dimension 
three [33]. Witten's technique generalizes to all higher dimensions, as long as the 
manifold is spin [SI E] • In higher dimensions, a complication arises in the Schoen- 
Yau argument due to possible singularities of minimal hypersurfaces. Two different 
approaches for handling this complication have been announced by J. Lohkamp in a 
2006 preprint [19] and by the last author in 2009. The first author generalized the 
spacetime E > theorem to dimensions less than 8 (without spin assumption) in his 
thesis (see [TT]). 

For earlier history of this problem, see [25]. The E > theorem is sometimes re- 
ferred to as the "positive mass theorem," but we prefer to refer to it more accurately 
as the "positive energy theorem," and reserve the phrase "positive mass theorem" 
for the E > \P\ theorem, which could also reasonably be called the "future time- 
like energy-momentum theorem." One can see that there is a natural relationship 
between these two theorems by "boosting" an initial data slice of a spacetime. In 
the case of vacuum initial data this approach can be made rigorous by using [8] to 
approximate the data by an initial data set which coincides with a Kerr data set out- 
side a compact set. The boost can then be done near infinity in the Kerr spacetime. 
In the last section of the paper we show that this argument can be generalized to the 
case in which the data satisfies the dominant energy condition. More precisely we 
show that a general initial data set satisfying the dominant energy condition can be 
approximatied by one which is identical with a Kerr slice near infinity and still sat- 
isfies the dominant energy condition. This gives an alternative proof of the positive 
mass theorem from the positive energy theorem. It is discussed in the final section 
of the paper. We emphasize that the proof presented in this paper is self-contained 
and gives a new proof of the positive energy theorem as well. 

We point out that our main theorem does not include a rigidity statement. The 
expected rigidity theorem, which has been proven in the spin case [Ml El Ej , states 
that if, in addition to the hypotheses in the main theorem, we have E = \P\, then 

= |P| = and {M,g) isometrically embeds into Minkowski space with k as its 
second fundamental form. 
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We briefly review the minimal hypersurface proof of the time-symmetric positive 
mass theorem in [251 122], since our proof is closely related. The argument is by 
induction on the dimension 3 < n < 8, and proceeds by contradiction. Suppose that 
there exists an asymptotically flat Riemannian manifold (M, g) with nonnegative 
scalar curvature -R^ > and negative mass E < 0. By a density argument [20], one 
may assume that (M, g) has harmonic asymptotics and positive scalar curvature. 
The harmonic asymptotics and E < assumptions imply that for large enough A, 
the coordinate planes = ±A are barriers for minimal hyp ersurf aces. Consider 
an (n — l)-dimensional vertical cylinder dCp of large radius p in the asymptotically 
flat coordinates. For each h G [—A, A], let Sp_/j be an area-minimizing hypersurface 
whose boundary is equal to the height h sphere on dCp. As long as n < 8, this area- 
minimizing hypersurface must be smooth. Each Sp must lie between the barrier 
planes = ±A, and the area is actually minimized over h by some 

hp e (—A, A). This means that this T^p^hp has the property that 



for any variation $t of ^p,hp that is equal to vertical translation along dCp. One 
can then extract a smooth subsequential limit of S^ /j^ as p — )■ oo. This Sqo is 
itself an {n — l)-dimensional asymptotically flat manifold with energy equal to zero. 
Moreover, T.^ is a stable minimal hypersurface of M. Not only does ([T]) imply that 
it is stable with respect to compact variations, but it is also stable with respect to 
variations that are (sufficiently close to) vertical translations outside a compact set. 
Using the (by now) well-established relationship between the stability of minimal 
hypersurfaces and scalar curvature, the stability of Sqo allows one to construct a 
conformal factor that changes the metric on Sqo to one with zero scalar curvature. 
The stability with respect to variations that are close to vertical translations is then 
used to show that the conformal factor must decrease the energy of Sqo, thereby 
violating the time-symmetric positive mass theorem in dimension n — 1. For the 
base case of the induction, when n = 3, one can show that the stability of Eoo 
combined with its asymptotics at infinity are incompatible with the Gauss-Bonnet 
Theorem. Note that when n = 3, choosing a special height hp turns out to be 
unnecessary. 

We now compare this with our approach to the spacetime positive mass theorem, 
which is essentially a generalization of the proof described above. In particular, it 
does not use the time-symmetric positive mass theorem as an input, as was done 
in the Jang equation approach of [27J. Again, the proof is by contradiction. Let 
(M",(yf, /c) be an asymptotically flat initial data set satisfying the dominant energy 
condition and E < \P\. By a density argument, which we give in Section [6| we may 
assume that (M, g, k) has harmonic asymptotics and satisfies the strict dominant 




(1) 
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energy condition /i > \J\. We may assume further that P points in the vertical 
direction —dn in the asymptotic coordinate chart. The harmonic asymptotics and 
E < \P\ assumptions imply that for large enough A, the coordinate planes = ±A 
are barriers for marginally outer trapped hypersurfaces (MOTS). Again, we consider 
an (n — l)-dimensional vertical cylinder dCp of large radius p in the asymptotically 
flat coordinates. For each h G [—A, A], the existence theory of [12] allows us to use the 
barriers to construct a MOTS Sp^/j whose boundary is equal to the height h sphere on 
dCp. As long as n < 8, this MOTS is smooth, because it has a A-minimizing property 
(see Section [2|. Moreover, Ep^/j lies between the planes = ±A and is stable in 
the sense of MOTS with boundary [16] (see Section |2|. Since MOTS are not known 
to arise from a variational principle in general, we cannot find a suitable height hp 
simply by minimizing some quantity. To this end, we introduce a new functional J-" 
on hypersurfaces with boundary on dCp, for which we can find hp G (—A, A) with 
(roughly) the property that 



dh 



^{^P,h) > 0. (2) 



h=ho 



Inequality (|2| in conjunction with the MOTS-stability of Ep /j then plays a role similar 
to that of pj) in the time-symmetric case. Note that the hp selected in the time- 
symmetric case by minimization would satisfy ^ in our more general argument. As 
before, we extract a smooth subsequential limit Eqo of Epjj^ as p — )■ oo. This Eqo is 
itself an {n — l)-dimensional asymptotically fiat manifold with energy equal to zero, 
and Eqo is a stable MOTS in M. Using the relationship between stability of MOTS 
and scalar curvature established in [17], one can construct a conformal factor that 
changes the metric on Eoo to one with zero scalar curvature. Finally, (|2| plays the 
role of ([T]) in establishing that the conformal factor must decrease the energy of Eqo, 
thereby violating the time- symmetric positive mass theorem in dimension n — 1. 

As in the time-symmetric case, the delicate height-picking argument is not required 
when n = 3. As before, we use the asymptotics of Eqo and MOTS-stability to produce 
a contradiction to the Gauss-Bonnet theorem. 



sets up the basic definitions and 
establishes the existence of the 



The structure of the paper is as follows. Section 2 
recalls some useful background material. Section 3^ 
MOTS needed for the proof. Section |4] completes the n = 3 case of the main theorem 
(Theorem [T|. The basic n = 3 argument, which is explained in detail in Sections 
3] and |4| was first sketched out in [23], Section 7.2]. Section |5] contains the parts of 
the proof that are specific to dimensions higher than 3, including the height-picking 
procedure alluded to above. Finally, the last section contains the relevant density 
theorem and its proof, which may be of independent interest. 
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2. Definitions, notation, and basic facts 

In this section we review some of the terminology used in this article, as well as 
some basic facts. 

Definition 1. Let S be a closed ball in centered at the origin, and consider the 
exterior region R" \ S. For any A; e N, p > 1, g e R, we define the weighted Sobolev 
space W'1'^{W\.B) as the space of Lebesgue measurable functions / on R" \ B such 
that 




Now suppose that M is a C*^ manifold such that there is a compact set i^T C M 
and a diffeomorphism M \ = R" \ S. Then one can easily define W'l'^{M) in 
a routine manner, by choosing an atlas for M that consists of the diffeomorphism 
M \. K = M."- \ B together with finitely many precompact charts, and then using 
the W_f(W' \ B) norm on the noncompact chart while using the W '"^ norm on 
the other ones. Note that the space W'lq{M) and its topology depend only on the 
diffeomorphism M^K = M" \ B. The definition of W^'^{M) also extends to tensors 
simply by considering their components with respect to these coordinate charts. We 
will usually write W^g for W'lg{M) when the context is clear. 

Definition 2. Consider the exterior region R"\S. For any A; e N, a G (0, 1), g G R, 
we define the weighted Holder space C!';(R" \ B) to be the space of functions / on 
R" \ S such that 

|a|<fe ^ \a\=k 

If M is a C'^'" manifold such that there is a compact set K and a diffeomorphism 
M \ ^ R" \ then one can define C^^(M) just as we defined W'lf{M) above. 

Definition 3. Let n > 3 and a G (0,1). An initial data set is a complete C'^'" 
manifold M" equipped with a Cf^" Riemannian metric g and a Cj^'^ symmetric (0, 2)- 
tensor k. On an initial data set, one can define the mass density fi and the current 
density J by 

J = divg k — d{tig k). 
We say that {M,g,k) satisfies the dominant energy condition if 

A* > \ J\9- 
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It is often convenient to consider the momentum tensor 

n = k - (tig k)g. 

Note that vr contains the same information as k since k = n — ^^^(tTg7!')g. 

Let p > n, q > and go > 0. We say that the initial data set {M,g,k) is 
asymptotically fla^ if there is a compact set C M and a 0^'°" diffeomorphism 
M \ K = \ B for some closed ball 5 C M", such that in these coordinates, for 
i,j = l,2,...,n, 

{g^J - 5.,, e ly'f (M" \ i?) X Wbl^iR^ \ B) 

and 

(/i, J) e C°'°_q^. 

We sometimes say that (M, fc) is asymptotically fiat of type {p, q, go, «)• 

Under these hypotheses, one can define the ADM energy E and ADM momentum 
P to be 

~ n 
P n 

Pi = 7 — TT lim / T^iii^ndHn'^, 2 = 1,2, ...,n 

where z/q = ^ is the Euclidean outward unit normal, the measure "Hq"^ is Euclidean 
Hausdorff measure, and Un-i is the volume of the standard unit sphere in M". 

Remark. Theorem [T] still holds if we allow multiple asymptotically fiat ends in the 
definition of initial data sets. We simply use the large celestial spheres in the "other" 
ends as barriers in the proof of Lemma |6| 

In order to execute our main argument, we require better asymptotic behavior. 

Definition 4. Let n > 3, and let [M^-^g, k) be an asymptotically fiat initial data set 
of type {p,q,qo,a). We say that {M,g,k) has harmonic asymptotics if there exists 
a C^'" diffeomorphism as in the definition of asymptotic flatness, as well as a €2'^^ 
function u and a Ca!"^ vector field Y such that for i, j = 1,2, . . . ,n, 

u{x) = l + a|xp-" + 02."(|a;|i-'^) 

Y,{x) = 6i|xp-'^ + 02'"(|x|^-") 

4 

2 

Tiij = M"-2 [{LYS)ij - (div5F)5ij] 
"'^Note that the term asymptotically flat has various mcompatible definitions in the hterature. 
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where a, bi are constants, and Ly is the Lie derivative. The 0^'°(|x|^ ") indicates 
that the function hes in the weighted Holder space C^!f„. 

Notation. Let {M"',g, k) be an initial data set, and let E"~^ be a C^'°' hypersurface 
with boundary in M", equipped with a global unit normal u. Let D denote the 
covariant derivative with respect to the ambient metric on M. (In contrast, we will 
usually use V when differentiating functions.) We define the second fundamental 
form By, and shape operator S^, of S using the convention 

By{X,Y) = {Sj:{X),Y) = {Dxi^,Y) 

for vector fields X, Y tangent to S, where the angle brackets denote the inner product 
g. We define the mean curvature scalar Hj^ to be the trace of 5*2. Using this 
convention, if u is the upward-pointing normal of the upper hemisphere in W^, the 
mean curvature is positive. We also define the expansion 

e+ = Hj^ + trs k, 

where tr^ k denotes the trace over the tangent space of E. If 9^ vanishes on all of E, 
then we say that E is a marginally outer trapped hypersurface, or MOTS for short |^ 
Note that the property of being a MOTS depends on the choice of normal. For any 
vector field X defined along E but not necessarily tangent to it, we define 

n-l 

divsX = 5^(De,X,e,), 

i=l 

where ei, . . . , e„_i is a local orthonormal frame of E. 

Notation. Given a vector field X defined along E, we can decompose X into its 
normal and tangential components 

X = ^v + X. 

Throughout this paper, whenever there is a vector field with the variable name X 
on a hypersurface E, the function ip and the vector field X are defined according to 
this relationship. We use rj to denote the outward unit normal of 9E in E. 

The expression for the linearization of the expansion stated in the following propo- 
sition generalizes the well-known variation formula for the mean curvature. See, for 
example, [T7] . 

Proposition 2. Let E""^ he a two-sided hypersurface with boundary in an initial data 
set {M"-,g, k), and let v be a continuous unit normal field along E. Let X G V{TM) 
be a vector field, and let be the flow generated by X . We can compute the 



While both the nomenclature and the acronym are slightly inappropriate in our context, we 
find it convenient to use these terms. 
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expansion of the push forward Ht ■= with respect to the unit normal that 

points in the direction of ^t*{^) one? pull it back to a function on S using $f. The 
derivative of this function in t at t = is denoted by D6^\y.{X) and one has that 

De+UX) = -A^if + 2(W^s, V<^) + (divE Wj, - \W^\^ + Q^)ip + v^e+, (3) 

where 

Qs = - /i - J{^) - l\kj: + 
Here, k-z denotes the restriction ofk to vectors tangent to S, and Wy. is the tangential 
vector field on S that is dual to the 1-form k{v, ■) along S. 

Notation. Throughout this paper, we will frequently drop the S subscripts when 
there is minimal ambiguity. Everything is computed with respect to the metric g 
unless stated otherwise, so we will not write g explicitly in any notation that we use. 
In particular, we use % to denote Hausdorff measure with respect to the metric g, 
unless indicated otherwise. 

Definition 5. Let S be a MOTS in an initial data set {M,g,k). We define the 
operator 

Lj:v := -A^v + 2{Wy, Vv) + (divs Wy: - \Wy\^ + Qy)v, (4) 

where f is a function on S. Although this operator is not self-adjoint, the Krein- 
Rutman Theorem shows that there exists a (Dirichlet) eigenvalue with minimal real 
part, which is called the principal (Dirichlet) eigenvalue of L^. Furthermore, this 
eigenvalue is real, and if S is connected, the corresponding eigenspace is a one- 
dimensional space generated by a C^'° principal eigenfunction that is positive on the 
interior of S [161 P- 3]. If the principal eigenvalue is nonnegative, we say that S is a 
stable MOTS. For a closed MOTS, this concept goes back to [IHl E]- It is easy to see 
that this generalizes the notion of stability of minimal hypersurfaces with boundary. 

Proposition 3. Let H be a stable MOTS in an initial data set {M,g,k). For any 
compactly supported function v onH that vanishes along dT,, we have that 

{\Vv\'' + QYV^)dU''-^>Q. (5) 

s 



This result follows from an argument in [17J, cf. the proof of Lemma 15 below. 

We state the following geometric variant of the strong maximum principle for 
ordered hypersurfaces that are subsolutions and supersolutions of the same prescribed 
mean curvature equation. We refer to [21^ Lemma 1], [3l Proposition 3.1], and [U 
Proposition 2.4] for similar results and partial proofs. It is important to pay attention 
to the choice of normal here, and a good example to have in mind is the following: 
Consider a sphere of radius 2 that is tangent to a sphere of radius 1 in Euclidean 
space. Either the smaller one is enclosed by the larger one or lies outside it. The 
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(obvious) conclusion from the lemma in this simple example is that the larger sphere 
cannot lie inside the smaller one. 

Proposition 4 (strong maximum principle). Let g be a complete Riemannian 
metric on M := B2~^{fd) x (—2,2) C M"", let F he a function on the unit sphere 
bundle S'^'^iTM), and let ui,U2 G C^iB'^-^Q)) be such that -1 < ui{x') < U2{x') < 
1 for all x' G -B"~-^(0,1). We assume that the hypersurfaces with boundary Ej := 
graph{ui) C M are such that Hj]^{x) < F(x, z/E^(a;)) for all x G Si and H^^i^x) > 
F{x, 2/22(3;)) for all x G S2. Here, the mean curvature is computed using the upward 
pointing unit normals. If Si and S2 intersect at an interior point or are tangent to 
each other at a boundary point, then they must be equal. 

Remark. We will use this proposition in the following context: We will have two 
open sets Qi C ^2 such that Si and S2 lie on dQi and 80,2, respectively, and 6^ < 
while e+^>0. 

Next, we recall a few facts about conformal changes of the metric. In general, 

4 

given two metrics on M related by g2 = u^^'^gi, the scalar curvatures are related by 
the formula 

R2 = u-'^^[r,u-'-^A,u). (6) 

-2 

Given any hypersurface S in M, its unit normals satisfy z/2 = M"-^ jy^^ while the mean 
curvatures are related by 

H2 = [h, + ^«-iV.,«) . (7) 

Also, if (M", ^fi) is an asymptotically fiat Riemannian manifold and u = l + a|xp~" + 
Q2,aQ^|i-n^ is C^'", then {M"',g2) is also asymptotically fiat, and the energies are 
related by the formula 

E2 = E^~ ^ hm / uV dU^' (8) 

= Ei + 2a. (9) 

Finally, the proof of the positive mass theorem in the time- symmetric case uses 
regularity and compactness properties of minimizing hypersurfaces. Since MOTS 
are not Icnown to obey a useful variational principle, we rely on the theory of almost 
minimizing currents: 

Definition 6 ([10]). Let {M"',g) be a complete Riemannian manifold, and let T be 
an integral fc-current in M. Let U C M be an open set such that spt((9T) fl [/ = 0. 
Then T is X-minimizing in U if, for every integral {k + l)-current X with support in 
f/, we have 

Mu{T) < Mu{T + dX) + \Mu{X). 
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Here, M[/ denotes the mass of a current in U. 

This almost minimizing property was introduced and studied systematically in 
PH] . In [12] it was shown that the A-minimizing property is a natural feature of 
the MOTS that arise in the existence theory developed there, despite the absence 
of a useful variational characterization. The elementary properties of A-minimizing 
boundaries that are needed in this paper are also summarized in [121 Appendix A] 
with further precise references to the geometric measure theory literature. 



3. Construction of MOTS 

Our proof of Theorem [T] will be by induction on dimension and contradiction. Let 
3 < n < 8, and suppose there exists an asymptotically flat initial data set (M", g, k) 
of type {p,q,qQ,a) satisfying the dominant energy condition, but E < \P\. For the 
case = 3, we will obtain a contradiction to the Gauss-Bonnet Theorem in Section 
|4| and for 3 < n < 8, we will obtain a contradiction to the time-symmetric case of 
Theorem [l] in dimension n—1. By the density theorem (Theorem 18 ), we can assume 



without loss of generality, that [g, k) has harmonic asymptotics and satisfies the strict 
dominant energy condition jj, > \J\g. Specifically, we can choose asymptotically fiat 
coordinates (x^, . . . , x") on M \ =^ \ B for some closed ball B such that on 
\ 5, we have, for i,j = 1,2, ... ,n, 

4 

gij U"~^6ij 

2 

jCij = [{LY5)ij - {d\vsY)5ij\ 
for some u,Y E Cl'"^^ satisfying 

u{x) = l + a|x|2-" + 02'°(|x|i^") 
Yi{x) = 6i|xp^" + 02'°(|x|^-"). 

Without loss of generality, we assume that P = (0, . . . , 0, — |P|). 



3.1. Existence of horizontal barriers. 

Lemma 5. Let {M"',g,k) be as described above. Then, for sufficiently large A, we 
have that 6'|^„^^| > and ^|^„^_y^| < where the expansion is computed with respect 
to the upward pointing unit normal. 

Proof. It follows from ^ that the coefficient of the function u is just a = ^. 

We claim that the IxP"" coefficient of K is h = — ^^Pj, i.e. that 6„ = ^^\P\ and 
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bi = for i < n. To see this, note that 



{n - l)ujn-iPi = hm / TTijulf (H-Lq ^ 

= hm f (r,,, + r,- , - (div5 Y)6ij) vi du^-^ 

= hm / (2-n)|x|i-"[6i(z/o),+&,M.-&fci^o% + 0(|x|-i)]z/^rf?{o"-i 
= -{n - 2)hiUJn-i- 

We claim that 

^{tn^A} = i/{."=A} + tr|,n=A|(fc) = (n - 1)(|P| - E)K\x\-^ + 0(|x|-"). (10) 
To see this, we use harmonic asymptotics and formula ([T]) to compute 

TT 2(n-l) -^-Ir) 

= ^^^^"'[(S - n)a|x|-"x" + 0(|a;|-")] 
= -2(n - l)a|x|-"a;" + 0(|x|-") 

To compute tr|^n=A}(^), first note that 

kij = [{LY6),j - ^(div^y)^^^] . 



So we have 



n-l 



ra-1 



71-1 

^ ui^.^^^ {-^^Y^^J,, + 0(|x|-")] since = 0\\x\'-^) for z < 

ij=l 

-y„,„ + o(|x|-") 

(n-2)6„|a;|-"x" + 0(|a;|-") 
(n- l)|P||x|-"A + 0(|x|-"), 



n 



completing the proof of (10). Note that (10) shows that that for large enough A one 
has 6'|^„^^| > 0. The proof that < is similar. □ 




3.2. Existence of MOTS. 

Notation. We now fix A large enough so that Lemma |5] apphes. For large p, we 
define Cp:= K\J . . . , x") | {x^Y + . . . + < p^} to be the region 

horizontally bounded by a vertical coordinate cylinder of radius p, and we define Cp^h 
to be the part of Cp lying between the planes = ±/i. Define T '■= dCpD{x"' = h}. 

Lemma 6. Let {M"',g,k) and A be as described above. For every sufficiently large 
p and all h e [—A, A] there exists a C^'" stable MOTS J2p^h in Cp^A whose boundary 
equals Tp^h o-nd which meets dCp transversely. Each Tjp^h is a X-minimizing (relative) 
boundary in Cp^2A where X depends only on \k\co. Moreover, we can choose {'^p,h}\h\<A 
so that if —A < hi < h2 < A, then Sp,hi lies strictly "below" ^p,h2 the sense of 
relative boundaries in Cp^\. 

Remark. The regularity (and later, convergence) of the MOTS Sp^/j is the only place 
in the proof of Theorem [T] where the assumption n < 8 is used. For n > 8, the lemma 
still holds except that the A-minimizing boundaries Sp are only regular away from 
a thin singular set. 

Proof. First observe that by the decay conditions on g and k, the coordinate cylinder 
dCp has 6^ > with respect to the outward normal and 6^ < with respect to the 
inward normal. 

We would like to solve the MOTS Plateau problem with boundary equal to Fp /i for 
each h G [—A, A]. Note that Tp^h divides dCp^2h into a top piece and a bottom piece. 
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According to [13] the MOTS Plateau problem is solvable if the top piece has > 
with respect to the outward normal of dCp^2k and the bottom piece has 6*+ < with 
respect to the inward normal of dCp_2K- By the observation above and Lemma [5| 
the two pieces of dCp^2k satisfy the desired trapping conditions, with the exception 
of the corners where dCp intersects {x" = ±2A}. See Figure 1. Intuitively, the 
corners do not cause a problem because the (singular) distributional mean curvature 
there has a favorable sign. Therefore one could "round-off" the corners as in Figure 
1. Alternatively, we observe that the proof in [T2] can easily accommodate the 
corners: Just as in that proof, we use the trapping of the cylindrical and horizontal 
pieces of dCp^2k to construct barriers that have appropriate blow-up behavior for 
Jang's equation, and then we simply combine these barriers as in [131 Lemma 4.1]. 
Specifically, there exists a C^'" family of MOTS {Sp/j}|h|<A with dT,ph = Fp/^ such 
that each Ep/j is a A- minimizing relative boundary in Cp^2k-, where A = A(|/c|ci)). 

Moreover, it can be seen from the construction in [12] that the regions Op^^ C Cp^2k 
bounded by the Sp^/j's can be taken to be ordered, so that ^p^hx C ^p^h-i whenever 
—A < /ii < /i2 < A. To see this, note that the supersolutions used in the 

construction of solutions to the regularized Jang's equation that lead to blow up in 
[T2| Lemma 4.1, bottom of p. 570] can be taken to be pointwise decreasing in the 
parameter so that the corresponding Perron solutions wfp/j are decreasing in h 
and hence their epigraphs are increasing. It has been shown in [H] that we may 
assume further that Sp^/i is stable in the sense of MOTS. 

Standard barriers for Sp can be constructed from the trapped boundary dCp by 
slight (C^)-inward perturbation above respectively below its boundary Fp /j, locally 
uniformly in (p, K)^ so that the angle (in the underlying Euclidean coordinate system) 
at which Hp meets dCp is bounded away from 0. For the special case of MOTS 
this inwards bending is explained in some detail in [121 Section 3]. Together with 
the A-minimizing property and Allard's boundary regularity theorem [1] it follows 
that near dCp the hypersurface Sp^/i can be written as a vertical C^'" graph above 
{x" = K). In particular, we see that Sp^/^ intersects Cp transversely. Higher regularity 
of the defining function — which solves a prescribed mean curvature equation — then 
follows from Schauder theory in a standard way. 

Finally, since each horizontal plane above = A has 0^ > and each horizontal 
plane below x" = —A has < by Lemma [s] it follows from Lemma |4] that each 
Sp /i actually lies in Cp^A- D 

Remark. We are grateful to Brian White for helping us clarify the following issue: 
When M is an n-dimensional C^'" manifold and is a complete Riemannian metric 
on it, then the standard interior AUard-type C^'" regularity of (almost) minimizing 
boundaries away from a set of Hausdorff dimension at most 8 holds. This was shown 
by J. Taylor in [30j (this part of the discussion in her paper applies to n-dimensional 
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manifolds). When the manifold is and the metric is Lipschitz, then this follows 
also from the work of R. Schoen and L. Simon [2l] (for almost minimizers this 
was pointed out by B. White in p. 498]). When the manifold is and the 
metric so that the Nash embedding theorem provides an isometric embedding of 
(M, g) into a high dimensional Euclidean space, then this also follows directly upon 
applying the Euclidean regularity theory as in [29]. In the preceding lemma, note 
that once we know that our surfaces are C^'", we can then apply Schauder theory to 
the MOTS equation to obtain C^'" regularity. The boundary regularity follows more 
easily because the metric is conformal to the Euclidean metric there. 

3.3. Convergence of MOTS. Although there is no reason to expect the family 
{Spjj}|/,|<A to be a C'^'" foliation (even when there are no topological obstructions), 
we can still prove a partial regularity result. Its proof is similar to the proof of 
regularity of the outermost MOTS (established in [3J for n = 3 and then in p!3] for 
3 < n < 7), but simpler because the two-sided A- minimizing property eliminates 
concerns about embeddedness. 

Lemma 7. Let {^p^h}\h\<A be as in Lemma^ For each ho E (—A, A], the upper 
envelope of {J:p^h}h<ho 

is a 0^'°^ MOTS with boundary V p^h^ which we denote by S^^^^. 
By convention we define Sp_A := ^p,-\- Moreover, limh/^ho ^p,h = 'E.ph 
topology. We define T^p^ho as the lower envelope of {T,p^h}h>ho for ho G [—A, A) and 
Sp^A := Sp.A o,nd note that analogous statements hold for these hyper surf aces. 

Proof. Fix ho G (A, A] . Let —A < hi ho as i ^ oo and pass the A-minimizing 
boundaries Ep /^^ to a subsequential limit S^ ^^^. Current convergence is automatic 
from the mass bounds, varifold convergence follows because there is no mass loss in 
limits of A-minimizing currents, and C^'" convergence follows from Allard's interior 
and boundary regularity theorems. Note that since the Sp^/j. are increasing, this 
limit is independent of the choice of subsequence, and hence it is really a limit of the 
original sequence hi. Moreover, the limit does not depend on choice of sequence hi 
for the same reason. □ 

Definition 7. We say that ho G [—A, A] is a jump height if T^p^ho does not equal 

By the previous lemma, T,p^h converges to Sp^/^g in C^'" as h ^ ho precisely when 
ho is not a jump height. It follows from the preceding lemma and Lemma |4] that 
ho is not a jump height if and only if the map h — )■ C^iVLp^h) is continuous at /iq, 
where £"(f2p^/i) is the volume of the enclosed region fip/^ defined in the proof of 
Lemma |6j Note that this implies that there are at most countably many jump times 
ho G [-A,A]. 

We also observe that the MOTS constructed in Lemma [6] may be used to construct 
complete MOTS with good asymptotics. 
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Lemma 8. For any choice of pj — )■ oo and hj G [—A, A], there exists a subsequence 
ofLp- h^. that converges in C^'" on compact subsets of M to a complete C'^'" properly 
embedded MOTS Soo- Moreover, there exists a constant c G [—A, A] such that outside 
a large compact subset of M , Hqo can be written as the Euclidean graph {x" = f{x')} 
of some C^'" function f{x') = c + 0^'"(|a;'p~") in the {x^, . . . ,x^~^,x"') = (x',a;") 
coordinate system. 

Proof. Existence of a subsequential limit Sqo and C^'" convergence follow as in the 
proof of Lemma [Tj Note that the limit Sqo is again A-minimizing. Since each T^p^hp 
lies between the horizontal planes x" = ±A, so does Eqo- The estimates below take 
part in the complement of a large ball B in M where we have harmonic asymptotics 
for the metric so that gij = U"^Sij where u = 1 + 0^'"(|a;p~"). It is not difficult to 
see from the corresponding property of Sp^ that the vertical projection of Sqo onto 
the plane {xn = 0} fl (M \ B) is surjective. Note that ([7| implies that the Euclidean 
mean curvature of is 0(|x|^~"). The A-minimizing property of Eoo gives rise to 
an explicit estimate of the form 0(|x'|~^) for the Euclidean area excess of in 
large Euclidean balls centered at points {x',x^) G and of radius \x'\/2. Together 
with the AUard regularity theorem (the version in [29| Theorem 24.2] is particularly 
convenient here), this estimate imphes that outside some large compact set, Sqo is 
the graph of a function f{x') such that |/(a;')| < A and f{x') = 0^''^{1) for some 

7e(0,i). 

Since Sqo is a MOTS we have that Hj]^ = —tT-£^{k). As in [221 P- 32], this 
translates into a prescribed (Euclidean) mean curvature equation for / via a confor- 
mal change ([t]). The initial estimate / = 0^''~'{1), together with a computation of 
tr^,^ k as in the proof of Lemma |5| shows that the Euclidean mean curvature of / is 
0^{\x'\^~^~'^) for some 7 G (0, 1). Standard asymptotic analysis as in [201 EI] shows 
that there exists a constant c G [—A, A] such that f{x') = c + 0^''^(|x'p~"). (Note 
that there is no logarithm term when n — 1 = 2 because / is bounded.) Repeating 
the above analysis with this information shows that f{x') = c + 0^'"(|a:;'|'^~"), as 
asserted. □ 

Corollary 9. When n > 3, the hypersurface G M in Lemma^is asymptotically 
flat and has zero energy with respect to its induced metric goo- 

The following lemma is a simple consequence of Proposition |3] and the fact that 
Soo is a limit of stable MOTS. 

Lemma 10. Let Sqo be a complete MOTS whose existence is established by Lemma 

~~ rl,2 



a\ For any v G WsL^ (^oo); we have 

[ (|Vi;p + Qs^t;2) >0. 
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We omit the proof because it is strictly simpler than that of Lemma 17 in Section 



5.3 



Specifically, the proof of Lemma 17 becomes a proof of Lemma 10 by simply 
replacing Z by an arbitrary compactly supported vector field on M and replacing 
the use of Lemma 15 by Proposition |3l 



4. The case n = 3 

We consider the base case n = 3 of our inductive proof. After the preparation of 
the previous section, the rest of the proof of the n = 3 case is essentially the same as 
for the time-symmetric case in [25], where minimal surfaces are replaced by MOTS. 

By Lemma Isl we can extract a subsequential limit Sqo of Sp o as p — )■ oo. Let YJ^ 



be the noncompact component of Sqo. Lemma 10 implies that 



I (^^vf ^QT.<y)dU^>^ (11) 
^ E^ 



for every v G ly^(Soo). Noting that Sqo has quadratic area growth, we can use the 

2 

"logarithmic cut-off trick" exactly as in [251 pa-ge 54] to approximate the constant 
function 1 on YJ^ by compactly supported functions in order to conclude that 



The strict dominant energy condition then implies that 

Kj:, dn' > 0, (12) 



where Kj]/^ denotes the Gauss curvature. On the other hand, just as in [22], the 
estimate {goo)ij{x') — 6ij = 0^{\x'\~^) implies that the geodesic curvature of d(T,'^ fl 
Cr) isn= i+0(r-2) while the length of (9(S'^na) is 27rr+0(l). The Gauss-Bonnet 
Theorem tells us that 

/ Kj,,^ dn" = 27rx(E'^ n a) - / t^dUK 



Combined with (12) and the asymptotics of d{Y'^ flCr) described above, for large r, 
we obtain 

< 27rx(S'^ n Cr) - 27r. 

Since S'^ \ is a graph for large r, we know that S'^ fl Cr is connected, yielding a 
contradiction. □ 
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5. The case 3 < n < 8 

Let 3 < n < 8. We assume, in addition to the indirect assumption in Section [3] that 
Theorem [T] holds in n — 1 dimensions. For the reasons described in the introduction, 
the proof here is substantially different from the proof in the time-symmetric case. 
The fundamental problem here is how to choose the MOTS Sp with the appropriate 
height h. As seen in Section |4| this is not an issue when n = 3. 



5.1. The functional J-'. In this section we consider a functional J-' that will be 
essential for our proof. In order to motivate the definition of J-", consider the time- 
symmetric case when k = so that the MOTS {T,p^h}\h\<A constructed in Section 



3.2 are minimal hypersurfaces. An important step in the proof of the Riemannian 
positive mass theorem when 3 < n < 8 [22] is to pick Sp,/ip with least area in this 
family. Suppose for a moment that the family {Ep^/i}|/i|<|A| is actually a C^'"" foliation 
of minimal hypersurfaces with a first-order deformation vector field X that is equal 
to dn at dCp. Then 



dh 



^"-'(Sm)= / (divs,,, X) 

= [ div2,,,(^z. + X)dH" 

(divs,., X) dn 
(X, 7]) dn 



p,h 

{dn,r^)dn^-\ 

If hp G (—A, A) minimizes areas as described above, then the first derivative in h 
of the integral above is nonnegative. This, along with the stability of ^p,hp among 
deformations that keep the boundary fixed, is all that is needed to finish the proof 
in the time-symmetric case p2] . 

We now return to the general case. Instead of using the area functional, which is 
not adapted for application to MOTS, we will build our proof around the functional 
described below. 
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Definition 8. Let S be a compact hypersurface in M whose boundary lies on some 
coordinate cylinder dCr- Then we define 



^(S)= / {dn,v)dH"-', (13) 

where r) is the outward unit normal of 9S in E. Note that, using harmonic asymp- 
totics, one can easily see that 

^(S)= / u-^Vodnr', (14) 

where ?7q is the n-th component of the unit normal r^o computed using the Euclidean 
metric, and Hq"^ denotes Euclidean Hausdorff measure. 

The barrier planes {x„ = ±A} give us a sign on J^(T,p^±\): 

Lemma 11. For any p sufficiently large, 

^(S^,_a) < < ^(S,,a). 

Proof. From Lemma [6] we know that Sp A hes below the plane {x"^ = A} in Cp. The 
strong maximum principle Lemma |4] implies that they cannot meet tangentially at 
their common boundary Fp^A- Hence the Euclidean outward unit normal of 9Sp^A in 



Ep^A satisfies viq > 0. The inequality J-'(Sp^A) > then follows from (14). The proof 



that J-'(Ep _a) < is analogous. □ 



Recall the definition of jump heights from Section 3.3 



Lemma 12. The function h i— t- J-'(Ep^/j) is continuous at every Hq G [—A, A] that is 
not a jump height. If Hq G [—A, A] is a jump height, then 

lim J^(Sp,,,) > J^(Sp,ftJ > lim J^{^p,h), 

where both limits exist, and at least one of the inequalities above is strict. In other 
words, there must be a downward jump discontinuity at every jump height. 

Proof. Let Hq G [—A, A]. Then by Lemma [7| Yimh/-hoJ^{^p,h) = J^i^Lp^ho) ^^'^ 
\imfi\ho J^(X'p,h) = J^iX'pM)- definition, if Hq is not a jump height, then both 
of these limits must equal J^{T,p^ho) 

Let Hq be a jump height. Since the family {T,p^h}\h\<A is ordered, it is clear that 
^p^ho li^s beneath T,p^ho, which lies beneath Sp^/j^. Since they all share the common 
boundary Tp^ho, we have (?7)q > t]q > {fj)Q, where ri,ri,f] are the outward normals of 
Tp^h^ in Ep f^^, Ep^ho, and Ep^ho, respectively. Since Hq is jump height, E^^^^ ^ '^pM^ 
so the strong maximum principle (Lemma H| implies that at least one of the above 



inequalities is strict. By the definition of J-" in (14) 
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where at least one of the inequalities is strict. □ 

Notation. Let Z be any vector field on M that is equal to 9„ outside some 
compact set. Given a hypersurface S equipped with a global unit normal z/, we use 
the decomposition 

Z = 0Z/ + Z, 

into normal and tangential components along S. 

We now compute the first variation of T . In view of Lemmas [TT] and [12] we may 
hope to find hp G (—A, A) such that the derivative oih^ T{Tjp^h) at hp (defined in 
a suitably weak sense) is nonnegative. 

Proposition 13. Let J] be a compact hypersurface in M whose boundary lies on 
some dCr, equipped with a global unit normal v. Let X be a vector field along S 
that is tangent to dCr at 812. As long as Z = dn along dC,., we have 

D^|e(X)=/ {cpVip + G{X),r,)dK'-\ (15) 

JdT. 

where 

G{X) = DxZ -DzX + {^H + divE X)Z - 05(X) - i^S{Z). (16) 

Remark. Of course, it does not matter what Z is away from dCy, and it is un- 
necessary to introduce Z at all except for the purposes of readability and for later 
application. 

Proof. Note that J^(S) = jQj^{Z, rj) d'H"'^'^. Let ei, . . . , 6^-2 be a local orthonormal 
frame for T(9S). We can differentiate Z, t], and the induced measure on 9S to 
obtain 

DJ^UX) = 

{DxZ, r]) + (z, {D^X, v)v - |^(/^e,X, r7)e,^ + (Z, r^) div^s X 

(The derivative of ?7 is computed by differentiating the orthogonality relations.) 
Along 9S, we introduce the decomposition Z = iprj + Z^ into components that 
are normal and tangential to 9S. The second term in the integrand of (17) is 

n-2 \ n-2 

Z, {D,X, u)u- 5^(De.X, r,)eA = (Z, u) {D,X, u) - ^^(Z, e,) {D^^X, r,) 
1=1 I 1=1 

= (P{Dr,i^u + X),u)-{DzaX,ri) 
= + {D,X,u)) - {DzaX,r]) 

= {cj^Vif, v) - {<pS{X),r]) - {D^aX, T]). (18) 



9S 



dU''-'. (17) 
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The third term in the integrand of (17) is 



(Z, r^) divas X={Z,7^) (div^ X - {D^X, r,) ) 

= {Z, 7]) i^H + divs X) - ^{D^X, 7]) 
= {{ipH + divE X)Z, T]) - {D^,X, rj). 



(19) 



Notice that the first term in the integrand of (jlTj), the first two terms of (18) and 
the first term of ( 19 ) combine to give 

(0V</^ + DxZ + {^H + divs X)Z - (pS{X),r]). (20) 



The remaining two terms, which are the last term of (18) and the last term of (19), 
combine to give 

-{DzaX,r]) - {D^r,X,r]) = -{D^X,r]) 

= -{D^{ipu + X),r,) 
= -{ipD^u + D^X,?]) 
= -{ipSiZ) + D^X,v). (21) 

□ 



The result follows from combining (20) and (21). 



5.2. Height picking and stability. The following lemma, whose proof we defer to 
Section 5.4, will stand in for the geometric inequality ([T]) that was available in the 
time-symmetric case. 

Lemma 14. Let S'^^ denote the component ofUp^h that contains the boundary Tp^h- 
For any large p, there exists hp G (—A, A) and a vector field X along S' . that 



is equal to Z = dn at (9E' 



^p,hp, such that ifi = {X, v) > 0, 

Do+y ^ (x) = o, 



and 



Dj^y , {X) > 0. 

-•3,0 



(22) 



(23) 



Consider h "path" of C hyp ersurf aces. // this path were differen- 

tiable in h (and if the Hp^h were connected), then the proposition would be easy to 
prove. By Lemma 11, we could find hp such that 

d 



dh 



> 



h=ho 



and choose X to be the first-order deformation of the family Ep at h 
the inequality above would become (23), and the fact that each Sp^/i is a MOTS 



would imply that L'6'+|s^,„p (X) 



hp. Then 



0. Unfortunately, T,p^h need not be differentiable 
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in h. In fact, if M has appropriate topology, the family Tip^h must contain jumps. 
From Lemma 12 we know that J-'(Sp/i) can only jump down at a jump height, so 
the presence of jumps does not cause any problem for finding hp as described above. 
Even in the absence of jumps, the lack of differentiability in h presents a technical 
challenge. 

Notation. For the remainder of this section, we will abbreviate S'^ by Sp. 



Lemma 14 allows us to conclude the following stability-like property, which the 
reader should compare to Proposition [3j 

Lemma 15. Let p he sufficiently large, and choose X and ip as in the statement of 
Lemma 14 Then for any any function v on Sp that is equal to (f) = (Z, v) at dUp, 
we have 

' {\Vv\'' + Qv^)dW^+ ! {G{X),ri)dn^-^ >0 (24) 



where 

G{X) = G{X) + (j)ipW. (25) 

Proof. We begin by following the argument in [T7]. Using equation (|3]), along with 
the positivity of (f, we have 

De+\^^{X) = -A^ + 2{W,V^) + {divW -\W\'^ + Q)^ 

= -Aip + (|V<^|2 + \W\'^ -\W - V<^n<^ + {divW - + Q)ip 

= -Aip + iVlogv^lV - \ W - Vlogv^lV + {divW + Q)ip 

= -{A\ogip)ip - \W - Vlogv^lV + (diviy + Q)ip 

= [div{W -Vlogip)](p - \ W - Vlogv?|V + Qv^- 

Now let V be any function on Sp that is equal to (p at 9Sp. We multiply the above 
equation by v'^(p~^ to obtain 

v^ip'^De+\j:^{X) = [diY{W -V\ogip)y -\W -V\og^\\^ + Qv'^ 

= dw{v'^{W - V logy?)) -{W - V\og!f,2vVv) 

- \ W -Vlog^l'^v^ + Qv'^ 
= diY{v\W -V\og^)) 

+ \{W -Vlog^)v\^ + |Vt;p - \{W -V\og^)v + Vt;|^ 

- Vlog(^|V + Qt;=^ 

= diY(v\W - Vlog(^)) + |Vt;P + Qv^ 

(26) 

-\{W -V\og^)v + Vv\^. 
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Together with equation ( 22 ) , this imphes 

< |Vi;p + Qv'^ + div{v'^{W - Vlog(^)). 
Noting that f = = <y9 at (9Sp, we can integrate this to obtain 

0< / {iVvl"" + Qv^)dW-^ + [ {v\W -V\og^),7])dK'-^ 



in-2 



< / {\Vv\^ + Qv^)dK'-^ + ! {G{X),r])dK'-\ 



where we used (15) and (23) in the third and fourth hues, respectively. 



□ 



5.3. Analysis of the complete MOTS limit. The strategy for the rest of the 
proof of Theorem [l] follows that of the time-symmetric case treated in [22] closely. 
By Lemma |8| there exists a sequence pj — )■ oo such that Sp, converges locally in C^'" 
to some Soo that has the properties described in Lemma [sj 

While Lemma [To] is strong enough to show that Sqo is conformal to a scalar-flat 
asymptotically flat manifold, it is not strong enough to control the change of mass 
under this conformal change. For that, we need to be able to allow test functions 
that are asymptotic to L In order to prove this, we will need to have uniform control 
over the Sp's as in the following lemma. 



Lemma 16. Recall the definitions of Z and from Section \5.1\ The following 
estimate holds uniformly in p. 

(|V0p + Qs>') dK'-'+ f (G(Z),r^)dH"'M =0(^-1). (27) 

Proof. We would first like to show that the estimates 

D9^\^^{Z) = 0{\x\~n (28) 

DJ^l^^ncAZ) = 0{r-') (29) 

hold uniformly in p. Using harmonic asjTiiptotics, formula ([T]) applies outside of a 
compact set to give us 



n-2 
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where Hq and z/q are the mean curvature and upward unit normal with respect to 
the Euchdean metric. Clearly, Hq and uq do not change under vertical translation, 
and Hq is bounded, therefore the decay of u implies that DH\Y.p{Z) = 0(1x1"""). The 
tr^p A; term is easier to handle, giving us (28). To derive (29), we simply use (14), 
the fact that uq does not change under vertical translation, and decay of u. 

To prove the lemma, we vary Sp in the direction Z, which is just vertical translation 
outside a compact set. We repeat the exact same argument from the first part of 
Lemma 15, except we use the region Sp \ Cr instead of Sp, the vector field Z instead 
of X, the function (p instead of if, and we take v = (p. Then on Sp \ Cr, equation 
(26) becomes 

4>D9+\j,^{Z) = dw{4>^W - 0V0) + |V0|' + Q(f)^ - \W\^(P^. 
Note that need not be positive in order for this to hold. Therefore 

(|V0|2 + Q02)rfH"-^+ / {G{Z),r])dn''-^ 

JSpxC,. J9(Sp\C^) 

+ DT\^^{Z) - DTl^^nc^Z) 



I Isp(^) + IW^IV) dW-^ 



where the last line follows from the (28) and (29), decay of W , and volume control 
coming from the A-minimizing property of Sp. □ 

The following lemma is the analytical consequence of our careful height-picking in 
Section 15.21 

Lemma 17. Let v he a function on Sqo such that f — 1 G Wl'^„ (Sqq). Then 

I {\Vv\^ + Q,,y) dn^-^>Q. (30) 

Proof. Following the notation of Lemma |8] we use coordinates x' on Sqo \ B where 
5 is a large compact subset of M. It is easy to see from Lemma [8] that 

Qe^ = |i?E^ J{y^J - \\k^^ + B^S = odx'l-"). 

Since the volume growth of is like that of IR"'^^, one can then see that 

/ {\yv\'' + \Q^Jv^)dH''-'<oo, (31) 
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whenever t> — 1 G Wl^{T,oo)- 

2 

We will write Tij for Tip. = T'p_^j^^ in this proof and use a j-subscript to denote 
quantities that depend on Sj, including j = oo. Whenever r < pj, we use rjj to denote 
the outward unit normal of fl Cr) inside Tj. (In what follows, we can restrict 
our attention to values of r such that dCr is transverse to each Sj.) Fix a vector 
field Z on M that agrees with dn outside a compact set. Let (p^o = (^^Sooi^) ^"^^ 
(f)j = (z/s,,^). Note that Lemma jsj implies that cj)^ - 1 = O^dx'p-") e ly^(Eoo) 

and that G'oo(Z) = 0(|x'|^"") where Goo is as (|25|. We have that 



(|V0oo|' + goo0L) 



n-l 



lim 

r— >oo 



n-l 



lim ( [ (I V0oor + Qoo0L) dn""-' + [ 

(|V0,f + Q,03 ci?^"-^ + 
lim ( / (|V0jf + Q,-0^) / {G 



9(SoonCr) 



n-2 



lim lim 



{G,{Z),r^,)d'H 



n-2 



a(SjnCr) 



,(Z),r/,)dH 



n-2 



>o, 



where the last inequality follows from Lemma 15, and the fourth equality follows 
from Lemma 16 Hence (30) holds for v = 0oo- Moreover, since the argument works 
for any Z that equals (9„ outside a compact set, (30) holds for all test functions that 
agree with (p^o outside a compact set. 

We now argue by density that the lemma holds for any function v such that 



G Vr^(Soo). Note that C^^'^{T^) is dense in Whi{J:^). Let u 



be a 



sequence of functions in C^'"(Soo) that converges to v ~ (poo in W^L^ (Sqq). It is now 
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straightforward to check that 

/n-l 



< hminf [ (\Vvi\^ + QooV^) dK 

i—^oo 



l—^OO 



hminf / [{\Vv\^ + Q^v^) + {2Vv ■ V{vi - v) + 2Q^v{v, - v)) 



where the cross terms vanish because of (31). FinaUy, we note that 0oo — 1 G Wl'^„ 

2 



to obtain the desired result. □ 

We are now ready to derive a contradiction to our induction hypothesis. Without 
loss of generality, we throw away all compact components of Eqo , noting that Lemmas 



10 and 17 still hold on the remaining asymptotically flat component. The strict 
dominant energy condition implies that Q^o < ^Rs^, so Lemma 10 implies that for 
any nonzero v G Wl'^ (^oo); we have 

Using the fact that n > 3, we have 

/ (iV^I' + > 0' (32) 

for all nonzero v G Wl'^ (Sqq). This inequality implies that the conformal Laplacian 

^ ' 2 2 



is an isomorphism. In particular, we can find some nonzero v G Wl'Y j^oo) such that 



Setting w = 1 + t>, we have 



- ^^Rj^^w = 0. 



By elliptic regularity, w G Cfj°, and by [201 Theorem 2], we must have w{x') = 
'.. + 0^'"(|x'|^~"). In particular, if g^o is the induced metric on Sqo, then by applying 
(6) in dimension n — 1, the conformal metric w^^Qoo on SJ^^ is an asymptotically 
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flat Riemannian manifold of zero scalar curvature. Using ([8j) in dimension n — 1, the 
new zero scalar curvature metric on has energy 



= 0-^^f [ (iVwl'^ + wAw)d'H 



(n-3)w„_2 

rl,2 



Now since w — 1 E W-iL„ (Sqq) and is not zero, we can use Lemma 



n-l 

o; V i^ w a n 

plus the short 



argument used to derive ( 32 ) to see that 



17 



and hence E{w"-^ Qoo) < 0. This contradicts the time-symmetric case of Theorem [l] 
in dimension n — 1, which is our induction hypothesis. 

5.4. Proof of Lemma 1141 We define 

hp = mf{h e [-A, A] I J^{^p,h) > 0}. (33) 

Our goal for this section is to prove that this choice of hp satisfies the conclusion of 
Lemma 14 For now let us assume that Tip^hp is connected. 

Using Lemmas 11 and 12 it is easy to see that hp exists, lies in (—A, A), is not 
a jump height, and that Tlj^p^hp) = 0. Although the path h — )■ Sp /i of C^'" hyper- 
surfaces need not be differentiable in h at hp, the continuity of the map h i— )■ Ep /j 
at h = hp allows us to use the inverse function theorem to describe the family Ep /j 
precisely for h near hp, as we will see below. We note that the work of B. White 
[31] on the moduli space of minimal submanifolds with boundary in M" is useful for 
understanding the family Ep^/^. However, since we do not need the full power of [31], 
we choose to use a simpler approach similar to that of [15] in our Case 2 below. 

Since T,p^hp is transverse to dCp, we can find a C^'° vector field dr such that 
(dr,!^) > along Ep/i^, and dr = Z = 9„ at Tp^^^. By integrating dr, there exists 
a relatively open neighborhood U of Ep/^^ in Cp and a C^'° diffeomorphism F : 
'^P,hp X {—6,6) — > U such that F{-,0) is the identity map, and F^:(-^) = dr- In 
particular, F{-,t) maps dT^ph = Lp/j to T{hp + r) by vertical translation. 

Any C^'"" function w : Ep^/j^ — > {~6,6) defines a graph in Ep/^^ x (—6,6), whose 
image under the diffeomorphism F is a C^'° hypersurface in U which we denote 
graph[w], whose boundary lies on dCp. Since hp is not a jump height. Lemma [T] tells 
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US that for each h sufficiently close to hp, each T^p^h = graph [w] for some C^'" function 
w. Moreover, the function w converges to in C^'" as h approaches hp. 

Recall that the operator L^.^^^ on Sp^^^ defined by (|4) has a principal Dirichlet 
eigenvalue, which is nonnegative since Sp /j^ is stable. We consider two cases, de- 
pending on whether this principal eigenvalue is positive or zero. 

Case 1: The principal eigenvalue of L^.^^^ is positive. 

Define the map^^ : Co''°(Sp,;,J x M CO'"(Sp,;,J x M by 

Using (g and Q, we see that 

Since Ly,^^^ '■ Cq'" — )■ C"'" is an isomorphism, so is D'^\(Qfly The inverse function 

theorem tells us that there exists a function $ : (— e,e) — )■ Co'"(Sp^/ip) such that 
is the unique small solution of the equation = (0,s). That is, it is 

the unique C^'"-nearby MOTS whose boundary is V p^^^^g- So by continuity of Sp/j 
at h = hp and the uniqueness, we must have that Sp /j^^^ = graph[$(s) + s\ for small 
s. Therefore h ^ '^p,h is actually a map into the space of C^'" hypersurfaces at 



h = hp, so that the straightforward argument described in Section 5.2 applies. Specif- 
ically, if we take X = (^ + 1)9,- G C^'°, then X is a ffist-order deformation vector field 
for the family Ep^^. Since each Ep^/j is a MOTS, D9+\y^,^^{X) = ^9^^Jh=hp = 0. By 
construction of hp, we have DJ'Iy^^^^X) = J^J^(Sp /i) |^_^ > 0. Finally, the normal 
component of X is nonnegative on Ep /j and positive on Fp by construction, and 
hence positive everywhere by the strong maximum principle applied to the linear 
equation L-^^,^ = 0. This completes the proof of Lemma 14 for Case 1. 



Case 2: The principal eigenvalue of -^s^^p ^^^^ 



In this case, since Sp is connected, L-^^ has a one-dimensional kernel in Cg'" that 
is generated by a function that is positive away from the boundary. The same is then 



true for the adjoint ^ . We define the map ^ : Co'"(Sp,;,J xM^ — ^ C°'°(Sp,ftJ x^^ 



by 

m{w, K, s) = (^-K + ^Japh[^+«]' V{graph[w + s]), s 



'Technically, it is only defined on some neighborhood of the origin. 
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where V(T,) denotes the "signed volume" of the region of Cp lying above Sp^^,^ and 
below S. Then 

L>^|(o,o,o)K,«:',s') = lL^^,(iw' + s')dr,i^)-K', [ {{w' + s')dr,y)dK'-\s 

We claim that D\l'| (0,0,0) is injective. If D\E'| (0,0,0) ("W^'? 1^' 1 s') = (0, 0, 0), then obviously 
s' = 0, and since the image of Cq'" under L-^p,^^ is orthogonal to the kernel of 
^ , k' = 0. Finally, {w'dr, v) is in the kernel of L-£^^^ and has zero integral, and 
is therefore zero, proving the claim. Since D^'l (0,0,0) has index zero, it is also an 
isomorphism. So the inverse function theorem tells us that there exists a function 

($1, $2) : (e, e) x (e, e) Co''"(S,,,J x M, 

such that s), $2(^,-5)) is the unique small solution of the equation 

vi>($i(e,s),$2(e,s),s) = (o,e,s). 

In particular, graph[$i(^, s) + s] is the unique C^'"-nearby constant 6^ hypersurface 
whose boundary is Tp^hp+s and whose "signed volume" is ^. Define ^(s) := V{T,p^hp+s)- 
By continuity of Sp,/i a.t h = hp and the uniqueness, it must be the case that T,p^hp+s = 
graph[$i(^(s), s) + s] for small s. Note that since Tip^h^+g lies strictly above ^p,hp, 
the function + s must be positive. The complication here is that we do 

not know that ^(s) depends on s in any nice way. Below, we will see that this does 
not matter. 

By construction of hp, we can choose a sequence \ such that 

^(Sp,,^+,J > 0. 

Let us pass to a subsequence such that the unit vector in the direction {^{sk),Sk) 
converges. We consider two subcases. 

Case 2a: The ratio converges to some finite number as — i- 00. 

The basic idea here is that even though we cannot take derivatives as in Case 1, 
we can take subsequential limits of difference quotients instead. The hypothesis of 
this subcase implies that s^) + Sk] converges to some C^'" function 

Therefore we can take the limit of —9^ ^ ^_L 1 = to obtain 

Sk graph[<I>i(^(sfe),Sfc)+Sfc] 



We claim that the conclusion of Lemma 14 holds for X = ipd-r G C^'". The previous 
equation tells us that D6^\Y,pi^ (^) — 0' ^^"^ '^^ "^^^ see that the normal component 
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if) = f'{dr,i^) is nonnegative and equal to {Z, z/) > at Tp^hp- By the strong maximum 
principle for -^^Sph^, f is strictly positive. Finally, since /i^+s^) > 0, 

DJ^\^^ {X) = lim -[T{gTaph[<^Msk),Sk) + Sk]) - n^p,hp)] 

' fc— !>oo 

= — [J^(Sp,/, +,J - )] 

fc— !>oo Sfc 

> 0. 

Case 2b: The ratio converges to zero as A; — )■ oo. 

We will show that this is impossible. The assumption implies that ^j^[$i(^(sfc), s/c) + 
Sk] converges to some C^'" function that is zero at Tp^hp- As in Case 2a, we conclude 
that 

Let X = ifdr G C^'". Again, its normal component ip = ifidr, v) is nonnegative, but 
this time it has zero boundary values. Since we also have Lj^^i^^ip = 0, the strong 
maximum principle tells us that either ip is identically zero, or else < at ^p,hp- 
To see why ip cannot be zero, take the limit of the equation ^^^P''^p+''k^ = 1 to obtain 
L ipdl-V^'^ = 1. Just as in Case 2a, 

= hm ^[J-(graph[$i(e(s,), s^) + s^]) - H^pMp)] 
> 0. 

On the other hand, note that since X = along Tp /i^, only the second term in the 
identity (17) contributes to DJ^\j]^^ so that 



DT\^ {X)= [ {Z,{D,X, 



u)u)(m 

P,hp 



n-2 



(PVrjipdn''-^ 



■n 

P,hp 



< 



since V^v' < and > on T p^^. This contradiction shows that Case 2b is 
impossible, finishing the proof of Lemma [Ml when is connected. 

Now suppose T^p^hp is not connected, ancHet S'^ denote the component of Sp^/^ that 
contains the boundary. Observe that since hp is not a jump height, S'^^^ converges 
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to ^'p^hp as /i — )■ /ip, and also that = J^ijlp^h) for all h G [—A, A]. 

Using these two facts, we can see that the entire argument above can be executed 
with '^'p^h^ in place of Sp^/j. 

6. The density theorem 

Recall that ir = k — [tig k)g. In this section it will be convenient to formulate 
everything in terms of vr rather than k. With a slight abuse of terminology, we will 
refer to {M,g, it) as an initial data set. The purpose of this section is to prove the 
following theorem. 

Theorem 18. Let (M", g, it) be an asymptotically flat initial data set of type {p, q, go, «); 
as defined in Section^ Assume that the dominant energy condition jj, > \J\g holds. 
For every e > there exists asymptotically flat initial data {g, vf ) on M of the same 
type as {g,TT) and with harmonic asymptotics (all in the same chart), such that we 
have 

\\9 - gWwi'^ < ^ "-^^ h - ^\\wLf_^ < 

the strict dominant energy condition 

> \ j\g 

holds, and 

\E-E\<e and \P - P\ < e. 

Remark. If we assume appropriate higher regularity for (M, g, vr), then (M, g, vf) will 
have the same regularity. This follows from applying Schauder estimates throughout 
the proof. Also, our proof shows that vr) can be chosen to be of type {p,q,q'o,(y) 
for any given gg > go- 

The proof of the density theorem consists of two general constructions. In Section 



6.1 we deform ((7, vr) to an initial data set that satisfies a s trict dominant energy 



inequality /i > (1 + 7)1-^^19 for some constant 7 > 0. In Section 6.2 we apply a cut-off 
argument as in ^ Theorem 1] to perturb to harmonic asymptotics, arguing that we 
can preserve the strict dominant energy condition in the process. Observe that even 



if we do not ask for strictness in the statement of Theorem [18} our proof still requires 
us to first perturb to strict dominant energy as in Section OA 

We first include a general fact that ADM energy and ADM linear momentum are 
continuous on the space of asymptotically fiat initial data sets, as defined in Section 

m 

Proposition 19 (c/. Proposition 2.4]). Let {g,T^) and {g,7r) be asymptotically 
flat initial data sets of type {p,q,qo,a). Given e > there exists 6 > depending 
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only on e,n,p, g, go, 



then 



\ 2,p^^i,p , and II (/i, J) - (/i, J)||^;i,o.i such that if 

-q -q — 1 —n — qQ/2 



g — g\\y/-i,p < 5 and IItt — 7r||^i,p < 5, 



\E-E\<e and \P - P\ < e. 



The argument is now standard and goes back to [261 P- 50] (for E only) and [H p. 
198] in the case of vacuum data. We include the proof for the sake of completeness. 
They key point is that continuity holds as long as we assume appropriate decay of 
the constraints. 

Proof. By the definition of the energy E and the divergence theorem, we have that 
2(n - l)a;„_i£' = lim / y^{gij,i - gii,j)i^idl-i'^~^ 

J\x\=r ~1 



J\x\=r 

~l" / / X 9ij,ij 9ii,jj) 

J\x\>r ~7 



for any r sufficiently large, and similarly for E. Note that J2iji9ij,ij ~ 9ii,jj) = 2/i + 
, since it differs from the scalar curvature Rg by terms quadratic in dkgij 



0{\x\ 

and TTij. Since (/i — /i) G W^'^_ it follows that there exists tq large and depending 



onlyone,n,p,g,go, ||/i-/i||^o.i ||(^-5, 7r)||^^2.p^^i,p ^, and ||(^-5, 7r)||^^2,p^^^i,p 
such that for all r > tq 



-i-9o/2 



q-1 



-q-1 



|a;|>r 



i9ij ,ij 9ii,jj)\dx < {n l)ti;„_ie. 



For the difference of the boundary integrals note that 

^ ^ [9ij,i 9ii,j i9ij,i 9ii,j)] ^0 d^Q 



\x\=r 



< 



2$^ / mg-gUdU 

i,j,k •^\^\=''' 



n-l 




< 2'Hr\{\x\ = r}) ^ sup \dk{g - gUi 

i,j,k 1^1"'' 

< 2{L0n-ir'''')n%r''"^Cp\\g - gW^^,,). 



X 



(34) 



POSITIVE MASS THEOREM 32 

1— — 

where Cp is the constant governing the continuous embedding W_f_g C C_i^g for 
p > n in Euchdean space. By choosing 6 > \\g — g\\yy2,p sufficiently small (depending 

on r) we can ensure that (34) is less than {n — l)uJn-i€ so that in total \E — E\ < e, 
as desired. The proof for the linear momentum is similar. □ 

6.1. Perturbing to strict dominant energy condition. We define the constraint 
map 

^g,7r) = (-2/i, J) = {-Rg + \71\l - -l_(tr^ 7,)2^ div, tt) . 

It is also convenient to define the modified Lie derivative CgY = LyQ — (div^ Y)g of 
a vector field Y. 

The proof of the following lemma is essentially the same as the proof of the sur- 
jectivity of D<l>|(g^^) in [HI Proposition 3.1]. 

Lemma 20. Suppose {g — 5, vr) G W'^f x W^f_g. Define the linear map A : W^'^ x 
Wit, ^ by 

A{h,w)=D<t>\^g,^){h,w) - {S),\h]Ji), 
where Z^$|(g.7r) is the linearization of ^ at {g,7c). Then A is surjective. 
Notation. For better readability, we define 



n - 2' 

which will appear frequently in this section. 
For {u -l,Y)e wit let 

g = u''g and n = u^'^'^in + CgY). 

Define the operator 

T{u,Y) = {-2fMu',Ju'/'^), 

where {ft, J) are the mass and energy densities of {g^Ti). Explicitly, the components 
of T are given by 

-2~^u' = '-^u-'AgU -R,-^, {tig TT + tr, CgY)' 
+ {\n\l + 2iCgY),in''' + \CgY\l) 

Jju'/' = {diYg CgY + divg tt) j + (ti - l)^M-iM,fc(7r + CgY)'; ^^^^ 

- 2"~^".i*^9(^ + ^9^)' i = 1,2, . . . ,n. 
Here, all indices are raised with respect to g. 
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Lemma 21. Suppose that {g - 5,7r) G VT^f x VT^f.,. For any (/, V) G VT^'IL, t/iere 
exist {v,Z) G one? symmetric {0,2)-tensors h,w in so that 

DT\^,,o)iv,Z) + D^\^g,^){h,w) = if,V, + \h]j{). 

If in addition (/, V) G C^Jn-qo /^^ some go > anc? a > 0, anc? {g, vr) G Cf^° x C/^", 
i/ien we /iat;e (f , Z) G Cj!?^. 

Proof. The linearization of T at (1,0) is 

DT|(i,o)(t;, ^) = (^A,t; - 4Zfc;,7r^'' - ^ tr, vr div, Z, 



(36) 



where all covariant derivatives and raising of indices are with respect to g. For 
< a < (ra — 2) and p > n, DT|(io) : Wl'^ — )■ W^f-a is Fredholm with index 



zero, cf. |5j. Because the linear map A defined in Lemma 20 is surjective onto 



„ we can find C'^'" compactly supported symmetric tensor fields {(/i^, w^)}^]^ 
whose images A{hk,Wk) span a subspace that complements the image of DT|(io) 
in W%^. It follows that for every (/, V) G there exist {v, Z) G W'^l and 

{h.w) G span{(/ifc, Wfc)}^]^ so that 

Dr|(i,o)(t;,^) + I^$|{<;,.)(/i,«^) = (/,^, + ^/i;JO- (37) 

This completes the first part of the proof. 

Now suppose that (/, V) G C°'"_qo for some go > 0, and (5',vr) G Cf^^ x C/^". We 



view the system (37) as (ra + 1) Poisson equations in i;, Zi with a decaying nonho- 
mogeneous term. Note that the contribution D$|(g^7r)(^, w) to the nonhomogeneous 
term is a compactly supported Holder function. Using [20l Theorem 2], which treats 
decay of solutions of the Poisson equation, and combining it with weighted Schauder 
estimates and a bootstrapping argument, we can conclude that (f , Z) G C^'^^- 

□ 

Theorem 22. Let (M", ^f, vr) he an asymptotically flat initial data set of type {p, g, go, a), 
as defined in Section^ Assume that the dominant energy condition /i > \J\g holds. 
For any 6 > 0, there exists asymptotically fiat initial data {g, vr) of the same type 
such that 

and for some 7 > depending on 6, 

/i>(l + 7)|J|g. 
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Proof. Choose a C'^'" positive function / such that / = near infinity, 

and let {v, Z) G €2'^^ and {h,w) G Cf'" be a solution of the system 



whose existence is guaranteed by Lemma 21 (with V = 0). Our goal is to show that 
for sufficiently small t > 0, the formula 

g = [l + tvy{g + th) and tt = {1 + tvY^^in + tCgZ + tw), (38) 

gives the desired initial data {g, tt) in the statement of the theorem. Since we clearly 
have \\g — gWi^^.p < S and ||7r — vf < 6 for small enough t, it suffices to show 

that /i > (1 + 7)1 J|g for some 7 > that depends on t. 
In the following we denote u = 1 + tv and define 

<^i{l + tv,tZ,th,tw) = (-2/iM^ Ju"/2)_ 

By Taylor expansion, 

+tv,tZ,th,tw) 
= $1(1, 0, 0, 0) + D$i|(i,o,o,o)(^^, z, h,w)+n 
= (-2/i, J) + tDT\^,^o){v, Z) + tD^^g^^){h,w) + 7^ 
= {-2^i,J)+tU,\h^,Jk)+'R, (39) 

where there is a minor abuse of notation in the last line, and the remainder term 
IZ = 7l{x, t) has the form 

n{x,t)=t / I (i,o,o,o)+ri(^,z,h,«,) - ^"^"1 1(1,0,0,0)] {v , Z , h, w) dr. 

Jo 

Claim. We have that 

\n{x,t)\ < ct\i + |x|)2-2- = o{t^f). 

where C is a constant that does not depend on x or t. 

Proof of claim. Clearly, 7l{x,t) = 0{t^), so it suffices to work outside a large ball 
B that contains the support of {h,w). Note that when x E M \ B we have that 

(1,0 fl,o)+rt(v, z,h,w){'v , Z , h,w) = DT(ifi-j^rt{v,z){v, Z) and hence 

7^(x, t) =t [DT\{i^o)+rt(v,z) - ^^1(1,0)] {v, Z) dr. 
Jo 
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From (35) we see that the hnearization of T at {u, Y) is given by 

-DT|(„,y)(f, Z) 



^(tr, TT + tr, C,Y) tr, CgZ + 2(vr + CgYf{CgZ) 



kl, 



div,(£,Z), + {n- l)^{vu-')4n + CgY)'; + {n - lf-u-\,{CgZ)'; 



The claim follows from substituting (1, 0) + rt(f , Z) and (1, 0) for (m, Y\ using that 
(w, Z) G ^2!^^, and estimating each term in an obvious way. □ 



From (39) we have 

u'ji = fi + tf + 0{t'^f) and u'''^Ji = Ji + t\h'lJk + 0{ff). 
In particular, for sufficiently small t > 0, we have 

u'^pL > + -/ everywhere on M. (40) 

We claim that u'^\J\g < + ^ for t > sufficiently small. Choose /i to be a C^'" 
nonnegative function whose support is larger than that of h. Then 

=u-%g^^ -tg^'^hi + Oit'h)) 

Therefore 

= ^g^3 _ tg^ky^ + o{eh)){j, + tikp, + o(tV))(^, + -^fc + o{ef)) 

+ 0(tVl'/i + tVl/ + ^'/i + tV') 
= |j|J + o(tVl/ + m 

2 .r ,2/2 



(1^1. + ^) -'i\J\,-% + o[e\j\f + t^n 



where we choose t > to be sufficiently small in the last line, proving the claim. 
Observe that the computation above explains the motivation behind the definition 
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of A. The ^h''jJi term in A is chosen specifically so that the coefficient of t above, 
which is the first order change in \J\g under the deformation, vanishes. 

Now fix t > small enough so that u^Ji > jj + ^ and u^\J\g < \J\g + ^. Our 
assumptions imply that sup^j ^ < oo. For x G M such that \J\g{x) 7^ it follows 
that 

fi _ u'fi + ^ \J\g + tf/2 ^ 

\J\S UV\S \J\9+tf/^- \J\g + tf/4- ^ 

where 7""'^ := 1 + | sup^,/ Since /i > we conclude that /2 > (1 + 7)] J|g on M, 
as desired. Finally, note that 7 depends only on t and supj^^ □ 

6.2. Harmonic asymptotics. We first show that any asymptotically fiat initial 
data set can be slightly perturbed so that 

g = vr = u'/^CsY, 

outside a compact set, for some choice of u and Y, meanwhile prescribing any con- 
straints (/i, J) that are close to the original ones. After that, we will show that if the 
imposed constraints decay rapidly enough, then [g, vr) has harmonic asymptotics. 

The first result is a modification of [S], Theorem 1], which treated the case n = 3, 
/i = /i = 0, and J = J = 0. 

Lemma 23. Suppose {g — S,tt) G W^'g x W^f_g. There exist Co, 60 > so that given 
(/i, J) G W^2_q with — /i, J — J)\\]yO,p < 6 < 60, there is an initial data set {g,Tf) 
of the form 

g = vr = u'/^CsY, (41) 

outside a compact set, for some {u — 1,Y) G Wr^, such that the mass and current 
densities of {g,Tf) are (/i, J). Moreover, 

11^? ^ 5'llw2,p < Co5 llvf - vrll^^i.p < CqS. 

—q —1—q 

Proof. For A > 1 large define the cut-off initial data 

{9\)ij = X\9ij + (1 - Xa)%, tta = Xatt, 

where Xx{^) = x{^/^) ^^id x is a C'^'" cut-off function on that is 1 on {|a;| < 1} 
and on {|x| > 2}. Note that || {gx — g,nx — ir) ||^y2,p^,^i,p — )■ as A — t- 00. 

In the following, we usually suppress the subscript A when the context is clear. 
Define 

g = u'g n = u'/\if + CgY), 
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where Y) G Wl'^. Let (/t, J) be the mass and current densities of {g, it). Define 

the map T(^,^)(u,F) = (-2/i, J) from + x Wl'^ ^ W^-f-g- The finearization 

of T^g^i,) at (1,0) is 

DT(^g,ir)\{i,o){v,Z) 

- AZk.fi^^ - ^ tr^ TT divg Z, 

t5 S S \ 

diyg{CgZ)j + {n- l)-v,kTT'; - -vj tr^ n - -(div^ 7i)jVj , 

where indices are raised and covariant derivatives are taken with respect to g. The 
map T(g^7r) : [Wl'^ +1) x Wl^ — j- W^f_g is defined analogously. Because g G (0, ra — 2) 
and p > n, DT(g^ii-)|(i^o) and DT(g_^) |(i^o) are Fredholm operators of index for A 
sufficiently large (see [5]). 

Let Ki be a complementing subspace for the kernel of i5T(g 7^)1(1 o) in W'^fx W^2_q. 
Since the linearization _D$|(g 7r) : W'^'g x VF^lf.g — )■ W^f-q is surjective (see the proof 
for n = 3 in [S], which can be generalized to n > 3) and because DT(g ,r)|(i q) is Fred- 
holm we can find C^'° compactly supported symmetric (0, 2)-tensors {(/i^, w^)}^]^ 
whose images D^\(^g^.„){hk,Wk) form a basis for a complementing subspace of the 
image of DT(c,7r)|(i,o) in W^2-q- Let K2 = spa.n{{hk,Wk)}k=i- We define the maps 
T(g,^),T(,,^) : ifi XK2 ^ W^°t, by 

T(g,^)(M, Y, h, w) = <l>{u'g + h, u'/^n + CgY) + w), 

and 

T(3,^)(M,r,/i,w) = ^{u' g + h,u''^{Tx + CgY) + w). 

Note that the maps T(g^i,),T(^g^.„) are continuously differentiable. Using that (^f, tt) 
converges to ((7, vr) in W'^^ x as A 00 it is easy to see that DT(^^^i,'^\(^u,Y,h,w) 

converges to -DT(g^7r) | {m.y,^,™) as A — )■ 00 locally uniformly in {u,Y,h,w) E Ki x K2 
in the strong operator topology. Observe that -DT(g ,r)| (10,0,0) is an isomorphism by 
construction. We conclude from the inverse function theorem that there exists So > 
such that for all A > 1 sufficiently large, restricts to a diffeomorphism defined 
on an open neighborhood of (1, 0, 0, 0) (independent of A > 1) in Ki x K2 and onto 
an open neighborhood containing the 14^°f_^-ball of radius 25o centered at {—2fi, J). 

Using that T(g,^)(l, 0, 0, 0) = (— 2yU, J) and that ||(/i, J) — (/t, J)||i4/0,p — )■ as A — > 

— 2 — q 

00 we see that ||(/i, J) — (/t, J)||i4/0,p < 26 provided that ||(/i, J) — (/i, t/)||^o,p < 6 
and A > 1 is sufficiently large. Hence if 5 G (0, 60) there exist {u — 1, y) G Wl'^ and 
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Qi,a compactly supported symmetric (0, 2)-tensors w) such that T(g^i,){u, Y, h, w) = 
(— 2/i, J) and such that 

for some constant Ci > that only depends on {g,7f). The lemma follows easily 
from this. 

□ 

The following proposition would be completely straightforward, but for a subtlety 
that arises in the case n = 3. 

Proposition 24. Suppose that {M'^,g,7i) is an asymptotically flat initial data set of 
type {p,q,qo,a) such that 

g = TT = u'/'^CsY, (42) 

outside a large ball B for some {u — 1,Y) G Wl'^. Furthermore, assume that qo > 1 
(rather than simply qo > 0). Then vr) has harmonic asymptotics in the sense of 
Section [B 



Proof. By (35), outside a compact set, {u — 1,1^) satisfies 

^--^A.^ - ^ (tr, CsYf + \CsY\ 



n-2 



-U 



AsY, + (n - l)-u-\k{CsY)'; - -u-\,tis{CsY) = Jju'/\ 

for J = 1, 2, . . . , 77,. By Sobolev embedding, {u — = 0^'"{\x\~'^). Just as in the 
proof of Lemma [21} we may view this system as + 1 Poisson equations for the 

|max(-2g-2,-n-go))_ XhcU USlug the fact that 



functions m, 1^, so that A.s{u, Y) = 0°'{\x\ 



q > Theorem 2 of ^Oj and weighted Schauder estimates imply that 



u 
Y 



l + a|a;p-" + 02'"(|xp-"-T) 



|2— n— 7N 



(43) 



6i|xp-" + 02'°(|x| 

for constants a, hi, for any < 7 < min(2g + 2 — n,qo), outside B. In particular 
{u — 1,Y) G C2^n outside B. For n > 3, repeating the above argument using 
(m— 1, Y) = 0^''^(|xp~"') yields the desired result, since max(2— 2n, — n— go) < —n—1. 
For n = 3 we argue as follows. 

We expand the nonlinear "source" terms above. We find that 

'2 - \LY9\^-2tTsiLgY){divsY) + 3idivsY)^ 



\CsY\ 



.«j=i 



2 



Y Y 



(div.r)^ 
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while 
Since 

we have 
and 



Yi 



\x\ 



(div5 Yf 



(6 ■ x) 



\X\ 



+ 



Thus the constraint equation for [i becomes 

2|6|2 3(6 



A5M 



2|x|6 



+ 0"(|x| 



as long as 7 < go- Now observe that 
A, 



3!/ 'tij j 


_ 25ii 




|a;|4 


1 


2 







\x\ 



for z,j e {1,2,3}, 



For appropriate choice of constants Cq and Cjj, we have 



A, 



u - Co 



1 XiXj 



\X\ 



\X\ 



0"(|x| 



Now apply [20l Theorem 2] and weighted Schauder estimates to obtain the desired 
result for u. The argument for Y is analogous. □ 



6.3. Proof of Theorem |18[ Using Theorem 22 and Proposition 19 we may reduce 
to the case where > (1 + ■j)\J\g- For / > 1 define ^i{x) := ^(y) where 



1 when |x| < 1 
''^ when Ixl > 2 



\x\ 



Note that J) — (/i, ^)||cO.« and hence J) — (/i, ^)||i4/0,p tends to zero 



-10 



2-q 



as / — )■ 00. Also note that J) G so that it has improved decay. 

By Lemma 23 and Proposition 24 there exist initial data sets {gi,7^i) with harmonic 
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asymptotics whose mass and current densities are J), and such that \\{g — 91,71- 
TT;) ||^2,p^^^i,p — > as / — !■ 00. In particular, H^f — gi\\co — so that 



-l-q 



ml=gr^fJ^J,=^f\J\l{l + oil)) as I -^00, 



and hence 



6^ > 6(1 + 7)1^1, > (1 + I) \^iJ\ 



an 



for / sufficiently large. This means that the strict dominant energy condition holds 
for {gi,Tri). The convergence of mass and linear momentum as / — )■ 00 follows from 
Proposition [191 □ 



Remark. One subtlety in the proof of Theorem 18 above is that the decay conditions 
for (/i, J) improve in the step where we are reducing to harmonic asymptotics. In 
fact, if we are content to preserve the dominant energy condition without requiring 
strictness, we can arrange for (/i, J) to vanish outside a large compact set by choosing 
^ with compact support. 

Remark. Finally we remark that it is possible to prove the positive mass theorem 
from the positive energy theorem. We argue by contradiction. Assuming the dom- 
inant energy condition and E < \P\, we can use the previous remark and [8j to 
approximate our initial data set by one which is identical with a suitably chosen 
slice in a Kerr spacetime outside a compact set, such that we still have the dominant 
energy condition and E < \P\. This makes it possible to further deform the data in 
the Kerr region to a "standard" slice with negative ADM energy, contradicting the 
positive energy theorem. Note that the positive energy theorem was proven in [27] 
for n = 3 and the proof generalized to the cases 4 < n < 7 in [11] . 
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